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ABSTRACT 

We present broad-band observations of the afterglow and environment of the short GRB 111020A. An ex- 
tensive X-ray Hght curve from Swift/XRT, XMM-Newton and Chandra, spanning ^100 seconds to 10 days 
after the burst, reveals a significant break at w2 days with pre- and post-break decline rates of ax.i ~ -0.78 
and ax ^2 ^ -1-7, respectively. Interpreted as a jet break, we infer a collimated outflow with an opening angle 
of 9j « 3-8°. The resulting beaming-corrected 7-ray (10- 1000 keV band) and blastwave kinetic energies are 
(2-3) X lO'*'^ erg and (0.3-2) x 10"*^ erg, respectively, with the range depending on the unknown redshift of the 
burst. We report a radio afterglow limit of < 39 /iJy (3a) from EVLA observations which, along with our find- 
ing that < i^x, constrains the circumburst density to «o ~ 0.01 -0.1 cm"^. Optical observations provide an 
afterglow limit of / > 24.4 mag at 18 hours after the burst, and reveal a potential host galaxy with / « 24.3 mag. 
The sub-arcsecond localization from Chandra provides a precise offset of 0.80" ±0.1 1" (Ict) from this galaxy 
corresponding to an offset of 5 - 7 kpc for z = 0.5 - 1 .5 . We find a high excess neutral Hydrogen column density 
of (7.5 ±2.0) X 10^' cm"^ (z = 0). Our observations demonstrate that a growing fraction of short GRBs are 
collimated which may lead to a true event rate of > 100- 1000 Gpc""* yr"', in good agreement with the NS-NS 
merger rate of k, 200-3000 Gpc"^ yr"^ This consistency is promising for coincident short GRB -gravitational 
wave searches in the forthcoming era of Advanced LIGO/VIRGO. 

Subject headings: 



1. introduction 

Observations of the temporal and spectral evolution 
of short-duration gam ma-ray burst (GRB; Tgo < 2 sec; 
iKouveliotou et alJll993h afterglows are crucial to our under- 
standing of the basic properties of these events: their en- 
ergetics, parsec-scale environments, and geometries. From 
observations over the past 7 years, we now know that short 
GRBs have is otropic-equivalent energies of ~ 10^"- 10^^ erg 
( Bergeiil2007l) and circumburst densities of ^ 10"^- 1 cm~^ 
JSoderberg et al."2006'; 'Panaitescu"2006; 'Sti'atta et al.l 120071: 
[Perley et al. 2009b; Berger 2010; Fong et al. 201 ll); however 
these ranges are based on only a handful of events. The 
geometry, or degree of collimation, is the least constrained 
property, but is of particular interest because it directly af- 
fects the true energy scale and event rates. These param- 
eters aid our understanding of the explosion physics, the 
nature of the progenitors, and the potential detectability of 
short GRBs as gravitational wave sources. In particular, 
knowledge of the true energy scale may constrain the mech- 
anism of energy extraction from the central engine and the 
ejecta composition : vv annihilation powerin g a baryonic jet 
(iJaroszynskil 1 1 9931: iMochkovitch et all 1 19931) or magnetohy- 
drodyna mic (MHD) processes in a magnetically-dominate d 
outflow jBlandford&Znaiekl 119771; iRosswog et akl l2003h . 
Significant impro vement on th e short GRB observed rate of 
> 10 Gpc"^ yr"' (iNakar & Gr anot 2007) will have a critical 
impact on estimates for coincident short GRB-gravitational 
wave detections in the era of Advanced LIGO/VIRGO 
(lAbadie et al.ll2010l) . 
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The opening angles (Oj) of GRBs can be inferred from 
temporal breaks in the afterglow light curves ("jet breaks"), 
which occur at the time, tj, when the Lorentz factor of the 
outflow is r(f,) » \/9i\ a later break correspo nds to a wider 
opening angle (ISariet al.ll 19991: lRhoadsl[T999h . Jet breaks in 
the light curves of long-duration GRBs have led to an open- 
ing angle distribution with a range of ~ 2 - 20° and a me- 
dian of 7°, leading to beaming-corrected energies of = [\- 

COS(^ , )1£^.iso 

200r, 'Friedma n & BloomI 
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Racusin et al. 20 09|). For short GRBs, there is m ounting theo- 
retical (Eichleret al. '1989'; Naravan et al. "1992) and observa- 
tional (Fong et al. 2010; Berger 2010; Church et al. 201 1 ) evi- 
dence that the progenitors are NS-NS/NS-BH mergers and nu- 
merous simulations of post-merger black hole accretion have 
predicted collimated outflows with 9 , ^ 5-20° ([Popham et" 
Il99l lAlov et alJl200l iRosswogl llOOS: Rezzolla et al."201 
up to several tens of degrees (Ruffert & Janka 1999b| 
iPopham et aT]|1999l;lRezzolla et al.ll201 ll) . 

However, the detection of jet breaks in the afterglow light 
curves of short GRBs has proved to be challenging. They 
can in principle be measured from optical or radio observa- 
tions, but there are several caveats that make this particu- 
larly difficult for short GRBs. First, the brightness of optical 
and radio afterglows are sensitive to th e circumburst densities , 
which are typically low, ~ 10"^ cm"-^ dSoderberg et al.ll2006h . 
Indeed, of nearly 70 short bursts detected by Swift, only 2 
radio afterglo ws ha ve been detected over the past^7 years 
dBergeret aLi 120051: ISoderberg et al.l i2006t IChandra & Fraill 
1201 Ih . Similarly, only ~ 30% of Swift bursts have detected 
optical aft erglows, with a typical brigh tness at < 1 day of 
« 23 mag ('Bergerl l2010l: iFong et 31112011 ). making long-term 
temporal monitoring nearly impossible with ground-based fa- 
cilities. Second, in the optical band there can be significant 
contamination from the host galaxies, which are g enerally 
brighter than the afterglows at > 1 day (lBergeiil2010h . 
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On the other hand, the X-ray afterglow brightness is inde- 
pendent of the circumburst densit y (as long as the dens ity is 
> 10"^ cm"^ and hence Vc > vx; i Granot & Sarill2002l) . and 
host contamination is not an issue. In addition, the well- 
sampled SwiftlXKl light curves from ~ 1 min to ^ 1 day 
provide an unambiguous baseline against which we can mea- 
sure a subsequent break. Therefore, it is no surprise that the 
X-rays enabled the discovery of the first jet break in a short 
GRB. The X-ray afterglow light curve of G RB 051221A ex- 
hibited a br eak at »5 days, leading to 9j «7° (Soderberg et al.' 
1^06; Burro ws et al.li2006) . Similarly, Chandra observations 
of GRB 050724 A out to 22 days pl aced a meaningful lower 
limit of 9j>25° (iGrupe et al.ll2006l) . consistent with a spher- 
ical explosion. Temporal breaks on timescales of < few 
hours were observe d in the XRT light curves of GRBs 061201 
(IStratta etal.ll2007h and 090510^ (De Pasquale et al. 2010). 
If these are interpreted as jet breaks, they would lead to 
9j ^ 1°; however, they also match the timescale and behav- 
ior of early breaks in long GRBs, which are no t due to col- 
hmation (Nous ek et ani2006t IZhanget all 12006': Lian g et all 
12007). Finally, there is tentati ve evidence for b eaming in the 
light curves o f GRBs 050709 (iFox etal.ll2005h and 061210^ 
(lBergeiil2007h : however, these are based on spar sely - s ample d 
light curves without a definitive break (e.g., I Watson et al.l 
|2006'). No other jet breaks in the light curves of unambigu- 
ous short GRBs have been reported to date^ and the lack 
of jet breaks in Swift/XRT observations out to 1 - 2 days 
can provide only weak lo wer bounds of 9j >2-6° (revised 
from ICoward et al1l2012l with more realistic density values; 
see Section|42|l. 

Against this backdrop, we present here the discovery of a 
break in the X-ray light curve of GRB 1 1 1020Aat St ^2 days, 
best explained as a jet break. We also present contempora- 
neous radio and optical limits on the afterglow, allowing a 
characterization of the broad-band synchrotron spectrum and 
constraints on the energy and density. In addition, we report 
the discovery of a putative host galaxy. A comparison of our 
X-ray and optical data may require an appreciable amount of 
extinction and the highest intrinsic neutral Hydrogen column 
density for a short GRB to date. Our results have strong im- 
plications for the opening angle distribution, and therefore the 
observed short GRB rate and true energy release. 

Unless otherwise noted, all magnitudes in this paper are 
in the AB system and are corrected for Galactic extinction 
in the direc tion of the burst using E (B -V) = .432 mag 
dSchlegel e t al. 1998; Schlafly & Finkbelniilloni). We em- 
ploy a standard ACDM cosmology with = 0.27, Ha = 
0.73, and Ho = 11 km s"' Mpc"'. 

2. OBSERVATIONS OF GRB 1 1 1020A 

2.1. Swift Observations 

GRB 111 020 A was detected on 2011 October 20 at 
06:33:49.0 UT by the Burs t Alert Telescope (BAT) on- 
board the Swift satellite (Gehr els et al.l l2004: Sak amoto et al.l 
1201 Ibl) . BAT located the burst at a ground-calculated position 
of RA=19''08"06.9^ and Dec= -38°0r50.3" (J2000) with 
2.1' accuracy (90% containment; ISakamoto et al.l201 lab . The 

^ G RB 090510 also exhibits a post- jet break-like decay in its optical light 
curve JNicuesa Guelb enzu et al.'2013. 

Please note that Berger (2007) erroneously refers to GRB 061006. 

^ A jet break was reported in the light curve of GRB 090426A 
INicuesa Guelbenzu et al.l 2011); however the characteristics of its envi- 
ronment and prompt emission are m ore similar to those of long GRBs 
ILevesgue et alJ20ia:IXin et al.l201 iT) . 



X-ray Telescope (XRT) commenced observations of the lo- 
cation of the burst at St = 72.8 s (where St is the time after 
the BAT trigger) and detected a fading X-ray source (Sec- 
tion |22|i. The UV-Optical Telescope (UVOT) began obser- 
vations of the field at St = 79 s but no corresponding UV or 
optical source was found within the XRT position. The 3(7 
li mit in the white fil ter, which transmits over A = 1600-8000 
A jPoole et al."2008'). is 20.3 mag (not corrected for Galactic 
extinction;|Oates & Sakamoto 201 1). 

The gamma-ray emission consists of a single pulse with 
a duration of Tgo = 0.40 ± 0.09 s in the 15 - 350 keV band 
classifying GRB 11 1020 A as a short burst (ISakamoto et al.l 
1201 lah . The spectrum is best fit with a single power law with 
index 1.37 ± 0.26 and a fluence of = (6.5 ± 1.0) x 10"^ 
erg cm"^ (15- 150 keV). Spectral lag analysis is not con- 
clusive, and there is no clear evidence for extended emission 
(ISakamoto etalJl2011h . 

2.2. X-ray Observations 

The XRT promp tly locate d a fading, uncatalogu ed X- 
ray source ( Evans et ar.ll2007l [2009; Sakamoto e t al. ,201 Ibl) 
with a UVOT-enhanced position of RA=19''08'"12.53' and 
pec=-38°00'43.8" (J2000) and an uncertainty of 1.6" 
(lOsborne et al.ll20TTI) . XRT observations of the field contin- 
ued until the source faded below the detection threshold at 
St Ri 3.5 days. 

We also observed the field of GRB 1 1 1020A with the Euro- 
pean Photon Imaging Camera (EPIC-PN) on-board the X-ray 
Multi-Mirror Mission (XMM-Newton) starting at St = 0.65 
days. With 13.5 ks of on-source observations, we clearly de- 
tect a source in the energy range 0.5- 10 keV, consistent with 
the Swift-XRT position. In addition, we obtained two sets of 
20-ks observations with the Advanced CCD Imaging Spec- 
trometer (ACIS-S; 0.3 - lOkeV) on-board the Chandra X-ray 
Observatory at St = 2.9 and 10.1 days to refine the astrometry 
and monitor the light curve evolution. We detect the X-ray 
afterglow in the first Chandra observation but do not detect 
any source at the same location in the second epoch. 

2.2.1. Data Analysis and Spectral Fitting 

We analyze the XRT data using the latest version of the 
HEASOFT package (v.6.11) and relevant calibration files. 
We apply standard filtering and screening criteria, and gener- 
ate a count rat e light curve following the p rescrip tions from 
Margu tti etlH (1201 Ol) and iMsrgutti et"aD (l2012l) . Our re- 
binning scheme ensures a minimum signal-to-noise ratio of 
S/N = 4 for each temporal bin. We analyze the XMM data us- 
ing standard routines in the Scientific Analysis System (SAS) 
V. 1 1 . We remove the first 5 ks of data due to high background 
contamination, giving a total exposure time of 13.5 ks. We 
extract count rates from a 20" radius aperture and the back- 
ground is calculated using 20" radius source-free regions on 
the same chip. We use the CI AO data reduction package for 
the Chandra data. For the first epoch, we use a 2.5" radius 
source aperture centered on the Chandra position and a back- 
ground annulus with inner and outer radii of 10" and 35", 
respectively, giving a source detection significance of ^ 30a. 
For the second epoch, we extract 1 count in a 2.5" aperture at 
the location of the source, consistent with the average back- 
ground level calculated from source-free regions on the same 
chip. We therefore take the 3(7 background count rate as an 
upper limit on the X-ray afterglow. 

To extract a spectrum from the X-ray data, we fit each of 
the data sets with an absorbed power law model (tbabs x 
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ztbabs X pow within the XSPEC routine) characterized by a 
photon index, F, and intrinsic neutral hydrogen absorption 
column, A^H.int, in excess of the Galactic column density in the 
direction of the bu rst, Nhmw = 6.9 x 1 0^"cm~^ ( typi cal un - 
certainty of ~ 10%: lKalberla et alj2005l:IWakker et 10201 lb . 
using Cash statistics. For XRT, we utilize data in the time 
interval 0.08-60 ks where there is no evidence for spectral 
evolution. We find an average best-fitting (C-stat,y = 0.86 
for 188 d.o.f.) spectrum characterized by F = 2.2 ±0.5 and 
A^H.int = (1.0±0.3) X lO^^cm-z atz = 0(Table[I]l. Uncertainties 
correspond to the 90% confidence level. Our best-fit param- 
eters are fully consistent with the autom atic spectrum fit pro- 
duced by 'Mang ano & Sakaniotol (1201 Ih . The XMM data are 
best modeled with a power law characterized by F = 2.0 ± 0.4 
and A^H.int = (0.65 ±0.22) x lO^^ cm'^ (C-stat^ = 1.0 for 256 
d.o.f.), consistent with the XRT model parameters. We also fit 
the first epoch of Chandra data and the resulting parameters 
are consistent with those from XRT and XMM, albeit with 
large error bars due to low count statistics (Table [T). Since 
we find no evidence for spectral evolution in the XRT data, 
we perform a joint XRTh-XMM spectral analysis to obtain the 
best constraints on F and A^H.int. The resulting best-fit model 
has F = 2.0 ±0.3 and A^e.mt = (0.8 ±0.2) x lO^^ cm'^ (90% 
C.I.; C-staty = 0.94 for 446 d.o.f.). Although the redshift of the 
burst is unknown, we note that F remains unchanged within 
its 1(7 value for z < 3 and we find evidence for intrinsic A^e.int 
in excess of the Galactic value at 6.5(7 confidence. The best- 
fitting spectral parameters for each of the three data sets and 
the joint fit are summarized in Table [T] 

Applying these parameters to the individual XRT, XMM, 
and Chandra data sets, we calculate the count rate-to-flux 
conversion factors, and hence their absorbed and unabsorbed 
fluxes (Table |2]i. Applying these conversion factors results in 
the X-ray light curve shown in Figure [T] 

2.2.2. Dijferential Astrometry 

In the absence of the detection of an optical afterglow 
(Section [2.3b . we use our Chandra observations to refine the 
SwiftlXKl position (1 .6" uncertainty) to sub-arcsecond accu- 
racy. We perform differential astromet ry be tween our Chan- 
dra and GMOS observations (Section 12.3b . To achieve the 
maximum signal-to-noise ratio, we combine both epochs of 
GMOS /-band observations and use SExtractor^ to deter- 
mine the positions and centroid uncertainty of sources in the 
field. Performing an absolute astrometric tie to 2MASS using 
^70 common point sources, we find a resulting rms value of 
crGMOS-aMASS = 0.17" (0.12" in each coordinate). 

To refine the native Chandra astrometry and determine the 
location of the X-ray afterglow relative to the GMOS im- 
age, we perform differential astrometry. We use CI AO rou- 
tines mergeall to combine the two Chandra epochs and 
wavdetect to obtain positions and \a centroid uncertain- 
ties of X-ray sources in the field. We also use wavdetect 
to determine the Chandra position of the afterglow from the 
first epoch and find a Icr centroid statistical uncertainty crx.ac; = 
0.08". We calculate an astrometric tie based on two X-ray 
and optically bright common sources and find weighted mean 
offsets of (5RA= -0.27" ± 0.06" and JDec= +0.05" ± 0.05" 
giving a tie uncertainty of ctcxo-gmos = 0.08". There are 
three additional common, but fainter sources. An astromet- 
ric tie using all five sources gives weighted mean offsets and 
errors of 5RA= -0.29" ± 0.15" and (5Dec= +0.06" ± 0.16", 

* http ://sextractor. sourceforge. net/. 



fully consistent with our results from using the two bright 
sources alone. We therefore use the astrometric solution from 
the two bright sources only. Applying this solution, we ob- 
tain a ChandraX-my afterglow position of RA=19*'08™12.49' 
and Dec=-38°00'42.9" (denoted by the blue circle in Fig- 
ure |2]i with a total la uncertainty of 0.20", accounting for 
the 2MASS-GMOS astrometric tie, GMOS-Chandra tie, and 
afterglow positional uncertainty. This position is consistent 
with the XRT position and is offset by 0.94" from the XRT 
centroid. 

2.3. Optical Observations and Putative Host Galaxies 

We initiated /-band observations of GRB 1 1 1020A with the 
Gemini Multi-Object Spectrograph (GMOS) mounted on the 
Gemini-South 8-m telescope on 201 1 October 21.01 UT {St = 
17.7hours). The data were reduced using the gemini/gmos 
package in IRAF. In a stack of 9 x 180 s exposures in 0.74" 
seeing and photometric conditions (Figure |2), we do not de- 
tect any sources within the enhanced XRT error circle or co- 
incident with the Chandra position. However, the outskirts 
of the XRT position are partially contaminated by emission 
from a nearby / = 17.7 mag star (SI) and a fainter star (S2) 
with / = 22.7 mag (corrected for A,- = 0.73 mag; Figure|2]l. We 
detect two additional nearby sources: a faint galaxy (Gl) lo- 
cated 2.8" away from the center of the Chandra position and 
a brighter galaxy (G2) with a 6.5" offset (Figure|2]i. 

To search for a fading optical afterglow, we obtained a 
second, deeper set of /-band observations (11 x 180 s) with 
GMOS at 6t = 1.74 days in 0.67" seeing . Digit al image sub- 
traction using the ISIS software package ( Alardl l2000l) reveals 
no variation between the two epochs inside or near the X- 
ray afterglow error circles (Figure 0. To calculate the upper 
limit on the afterglow, we add several point sources of vary- 
ing magnitudes between / = 24 - 26 mag around the position 
in the first epoch using IRAF routine adds tar. We perform 
photometry in the residual image in 2" apertures using the 
standard published /-band zeropoint for GMOS-S and place a 
3(7 limit on the afterglow of / > 24.4 mag (/v ^ 0.63 fJ,}y). 
We also perform photometry in a 1 .8" aperture for Gl and a 
2.3" aperture for G2, giving magnitudes of / = 23.9 ±0.2 mag 
and / = 21.9 ±0.1, respectively (Table|3]l. 

In addition, we obtained r-band observations (3 x 360 s in 
0.62" seeing) with the Low Resolution Survey Spectrograph 3 
(LDSS3) mounted on the Magellan/Clay 6.5-m telescope con- 
current to the first epoch of GMOS observations {5t = 17.7 
hours). We easily detect G2, with r = 21.1±0.1 mag, but do 
not detect Gl to a 3a limit of r > 23.4 (corrected for = 0.99 
mag; TableO, where the zeropoint is determined from several 
standard stars at similar airmass. 

Since S 1 and S2 contaminate the Chandra position, we sub- 
tract their contribution using point-spread-function (PSF) sub- 
traction on the individual observations and a combined stack 
of the two GMOS epochs. We use standard PSF-fitting rou- 
tines in the IRAF daophot package. Modelling the PSF us- 
ing 4 bright, unsaturated stars in the field out to a radius of 3" 
(^ 40fwhm) from the center of each star, we subtract several 
stars in the field including SI and S2. The clean subtraction 
of these stars indicates a model PSF representative of the PSF 
of the field. We uncover a faint, mildly extended source (G3) 
on the outskirts of SI at coordinates RA=19*'08™12.43' and 
Dec=-38°00'43.07" (J2000). This source, which lies 0.80" 
from the center of the Chandra error circle, has a magnitude 
of / = 24.3 ±0.2 and is a potential host of GRB 1 1 1020A (Sec- 
tion l3.1l ). However, we cannot rule out the possibility that this 
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source is a faint star. 

2.4. Radio Observations and Possible Afterglow 

We observed the position of GRB 111020A with the Ex- 
panded Very Large Array' (EVLA) beginning on 1 October 
20.95 UT {5t = 16.1 hours; Program IOC- 145) at a mean 
frequency of 5.8 GHz with a total on-source integration time 
of 65 min. We observed 3C295 and J1937-1958 for band- 
pass/flux and gain calibration, respectively, and used standard 
proced ures in the A stronomical Image Processing System 
? AIPS : lGreisenll2003:) for data calibration and an alysis. With 
the n ew wideband capabilities of the EVLA dPerlev et al.l 
l20Tlh . our data have an effective bandwidth of ^ 1 .5 GHz after 
excising edge channels and data affected by radio frequency 
interference. The low declination of GRB 11 1020 A and the 
compact D configuration of the array at the time of observa- 
tion caused significant shadowing and required the removal of 
7 out of 27 antennas (the north arm of the EVLA). 

Taking into account the highly-elongated beam (33" x 7" 
with a position angle of 170°), we detect a 3.7a source 
with an integrated flux density of 48 ± 13/iJy located 
at RA=19''08™12.40^ Dec=-38°00'41.2" (dRA= 1.1", 
(SDec= 3.6", 1(7 uncertainty), consistent with the Chandra 
position. The position, peak flux and integrated flux of the 
source are consistent regardless of our choice of weighting, 
or if we analyze the upper and lower sidebands separately. 
However, despite the statistical significance of the detection, 
we cannot completely rule out sidelobe contribution from 
nearby bright sources in the field due to the low declination 
of the burst. Therefore, we conservatively adopt a 3(t upper 
Umit of 39 ^Jy on the radio afterglow of GRB 111020A for 
our analysis. We note that if the source is indeed real then 
upper limits inferred from the radio data can be treated as 
actual values. 



3. RESULTS 

3.1. Galaxy Probabilities of Chance Coincidence 

To assess which galaxy is the most proba ble host of 
GRB 1 1 102( )A we adopt the methodology of Bloom et alJ 
(120021) and iBerged (1201 Ol) to calculate the probabihty of 
chance coincidence P{< 5R) at a given angular separation 
SR. We determine the expected number density of galax- 
ies brighter than a measured apparent mag nitude, m, using 
the results of deep op tical galaxy surveys (.Hogg et al.lll997l: 
iBeckwith et ai]|2006l) : 



cr(< m) = 



1 



0.33 xln(lO) 



X 10' 



,{).33(m-24)-2.44 



(1) 



Then the p robability of chance coincidence is given by 
(iBloom et aL,200 2) 

P«SR)=l-s-''^^'^>''"^^"'\ (2) 

We calculate P{< SR) for each of the three host galaxy 
candidates (Figure O, and find that G3 is the most probable 
host of GRB 1 1 1020A with P(< SR) = 0.01, while for Gl and 
G2, the values are P(< SR) = 0.10 and 0.12, respectively. 



' Newly renamed the Kai'l G. Jansky Very Large Array. 



3.2. X-ray Light Curve Fitting and a Jet Break 

The temporal behavior of the X-ray afterglow flux is char- 
acterized by a steady power-law decline until St ^ 2 days 
when there is a significant steepening in the light curve (Fig- 
ure [Til. A single power law model with a decline rate deter- 
mined by the X-ray data at early times (f < 2 days) provides a 
poor fit to the late-time data (dotted line in FigurefTll; in partic- 
ular, it overestimates the Chandra detection and upper limit. 
To quantitatively assess the shape of the X-ray light curve, we 
therefore invoke a broken power law model, given by 



Fy = Fx 



(3) 



where Fx.o = 2^I^Fx(t = tj), ax.i and ax.i are the power law 
indices pre- and post-break, respectively, tj is the break time 
in seconds, and 5 is a dimensionless smoothness parameter 
that characterizes the sharpness of the break. We perform 
a three-parameter x^-grid search over Fx.o, cnx.i and tj. If 
we use a relatively sharp break (e.g. s -10), the Chan- 
dra 3(7 upper limit constrains ax.2 ^ -1.7. If we allow for 
a smoother break (e.g. s « -1), ax.i remains unchanged 
but the break occurs at later times (tj « 4 days) and ax,2 
is required to have a steeper value of < -2.2 to accommo- 
date the Chandra upper limit. This scenario generally pro- 
vides a poorer fit to the last Chandra and Swift/XRT points. 
We therefore adopt the sharp-break scenario. Fixing .? = -10 
and ax, 2 = -2.1, we find a best-fit broken power law model 
characterized by Fx(tj) = (1.36 ±0.45) x lO"'-' erg cm"^ s"', 
axA = -0.78 ±0.05, and tj = 2.0 ±0.5 days (la, xl = 1-1 
with 15 d.o.f.; Figure[T]i- This best-fit model is shown in Fig- 
ure [T] The best fit parameters are independent of our choice 
of ax.2 between -1.7 and -3. We also note the presence of 
a slight flux enhancement relative to the power law decay at 
St oi 3 X 10'* s (0.35 days). If we remove these points from 
our fits, the resulting best-fit parameters remain unaffected. 

The required change in the temporal index is Aan 
0.9. There are several possibilities that can explain breaks 
in GRB afterglow light curves. The first scenario is the 
transition of the cooling freq uency across th e band, but 
this only predicts Aa = 0.25 dSari et all 1 1998b . An alter- 
native possibility is the cessation of energy injection, ei- 
ther from refreshed shocks or a l ong-lasting central en- 
gine (e.g., Rees & Meszaros] Il998t ISari & M6szarosl 120001: 
IZhang & Meszaros 2002), which has been used to explain 
the termination of a shallow decay or plateau phase in the 
X-ray and optical light curves of several long GRBs. How- 
ever, these cases all exhibit earlier temporal breaks at ^ 
10^- 10"^ se c with Aan ^ 0.7 (ax.i ~ -0-5, ax.2 » -1.2; 
iNousek et a l. 2006; Zhang et al. 2006; Li ang et al.ll2007l) . At- 
tributing the break in GRB 1 1 1020A to the cessation of cen- 
tral engine activity would require sustained energy injection 
from the start of XRT monitoring to the break time, ^100 
seconds to 2 days, whereas the timescales of energy injec- 
tion for long GRBs are < few hou rs (Nouseketal. 200^ 
IZhang et alJl2006l; iLiang et all [20071; iRacusin et al.1 12009). 
Single episodes of energy injection have also b een observed 
in two short GRBs: 051221A and 050724A fBergeretal 



in two short GRBs: 051221A and 050724A (IBergeretal 
2005 1 ; ISoderberg et ai1l2006l: iBurrows et al.ll2006llGrupe et al 
2006)). The light curve of GRB 05 1221 A, which exhibit! 



20061)). The light curve of GRB 05 1221 A, which exhibits 
a power law decay with index ax.i = -1-1, a plateau, and 
a return to the same power law (A an = 0), is interpreted 
as a single period of energy injection (ISoderberg et al.ll2006l ; 
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iBurrows et al.l2006l) . A super-imposed flare on the light curve 
of GRB 050724A with a single underlying decay index of 
cnx.i = -0.98 is also p ossibly relate d to late-t ime reactivation 
of the central engine (iBerger et al . 2005; Gr upe et all l2006l 
Figure[T]i. Neither of these light curves resemble the behavior 
of GRB 111020A, where the change in slope is substantially 
greater 

Yet another possibility to explain the break is a sharp 
change in the external density. However, models for den- 
sity jumps in a uniform medium (" Nakar & Gr anot' 2007) pre- 
dict that the density would need to decrease by greater than 
a factor of ^ 10"^ to account for the observed Aa^ > 0.9 
steepening. More realistic density contrasts of ^ 10 pre- 
dict Actmax » 0.4 in o ptical and X-ray afterglow light curves 
(iNakar & Granotll2007l) . 

Finally, we consider that the observed steepening is a jet 
break, when the edge of a relativistically-beamed outflow be- 
comes visible to the observer and the jet spreads laterally 
(ISari et al.111991 lRhoadslU"999l) . This model is often adopted 
to expl ain Aai2 ~J_Jn_the light curves of long GRBs (e.g., 
iFrail et al. 2001: iBlooniet a l. 2003; Ra cusin et a l. 2009) and 
has been observed in one other short burst, GRB 05 1221 A 
(Agp ~ 0.9, Figure fn I Soderberg et al.ll2006l iBurrows et aP 
|2006|) . Given the similarity in Aau and the timescales of jet 
breaks in both short and long GRBs, we conclude that the ob- 
served steepening in the light curve of GRB 1 1 1020 A is best 
explained by a jet break at tj = 2.0 ± 0.5 days. 

3.3. Afterglow Properties 

We utilize our radio, optical and X-ray observations to 
constrain the explosion properties and ckcumburst envkon- 
ment of GRB111020A. In particular, we adopt the stan- 
dard synchrotr on model for GRB afterglows (ISari et al.lll998[ 
iGranot & Saril |2002) which provides a mapping from ob- 
servable properties to the isotropic-equivalent kinetic energy 
(^K.iso), circumburst density (no), and the fractions of post- 
shock energy in radiating electrons (ep) and magnetic fields 
(eg). We use data at the time of the radio and first optical ob- 
servations {6t =17.7 hours), as well as the decay indices from 
the full X-ray light curve. 

First, we constrain the electron power-law index p, using a 
combination of temporal and spectral information. From the 
X-ray light curve, we measure ax.i ^ -1 .7 (Section | 3.2[) . For 
p = -ax. 2, appropriate for a spreading jet (ISari et al.ll999l) . we 
can then constrain p > 1 .7. To further constrain p and investi- 
gate the location of the cooling frequency, Uc, we compare the 
values ax.i =-0.78±0.05 and /3x=-1.04±0.16(/3x= 1-r; 
la) to the closure relations for a relativistic blastwave in a 
constant density ISM-like medium for p > 2, a typical envi- 
ronment expected for a short GRB from a non-massive star 
progenitor (Sari et al. 1999; Granot & Sai'i 2002). If i^c > i^x 
then the independently-derived values for p from the tempo- 
ral and spectral indices are inconsistent: p = 2.0 ± 0.07 from 
ax,i , and /:> = 3 . 1 ± 0.32 from /3x (errors are la). 

However, if i^c < vx we obtain p= 1.7 ±0.07 from ax.\, 
(jOranot & Sari 2002) which is consistent with the p value in- 
ferred from ax. 2, but yields a divergent total integrated energy 
in electrons unless a break at high energies in the distribu- 
tion is invoked. Although a flat electro n distribution (p < 2) 
is possible and not uncommon ie^'Dai & ChengI HqoI 
lPanaitescir & Kumar 2001; Ra cusin et a l. 2009), the standard 
relations for 1 < /:> < 2 yield p = 0.84±0.25 from ax.i. This 
solution is not self-consistent, and would also require an un- 
usually sharp break of A/? > 1.2 in the electron distribution. 



Continuing with the assumptions that Vc < vx and p > 2, we 
obtain p = 2. 1 ± 0.32 from Px, which is marginally consistent 
with the value inferred from the temporal index. Put another 
way, a - 3/3/2 = 0.77 ±0.30, w hich satisfies the closure rela- 
tion for Vc < Vx (ISari et al.lll998 ). We therefore conclude that 
Vc < Vx- We note that the spectral index is generally more 
reliable in the determination of p because it remains constant 
over time and is not subject to alternative processes such as 
energy injection or flaring. In this case, the same Px was also 
independently determined from both the XMM and XRT data 
sets (Table [T]). Therefore for the rest of our calculations, we 
take a reasonable value of p = 2.1 as determined from fix- 

We next determine a set of constraints on «o and ZiK.iso based 
on the X-ray flux density, radio limit, and the condition that 
Vc < Vx- First, we use the X-ray afterglow emission as a proxy 
for Ek.Iso assuming the X-ray emission is from the forward 
shock. For Vc < vx at the time of our broad-band observations 
(St = 17.7 hours), we use Fx = 0.032 ^Jy at^'x =2.4 x lO''^ Hz 
(1 keV), and p = 2. 1 to obtain ( iGranot & Sarill2002h 

£K,i,so « 2.2 X 10^"(1 +z)->e;> '"e^» °244« erg (4) 

where c/l.28 is the luminosity distance in units of 10^*^ cm. 
Next, we use isK.iso to constrain hq. Using our 3a EVLA limit 
of /v.rad ^ 39 /iJy, we can determine an upper limit on no 
under the reasonable assumption that our observed radio band 
(v = 5.8 GHz) is above the self-absorption frequency (vg < 

1/2 

Vmd < Vm \ fi/.rad PC »n ) at the time of observations. For this 
scenario (IGranot & Sarill2002l) . 

«o < 1.7 X IQ-^E-l'l.^H +z)-''X.2&4^'4'' cm-3, (5) 

where -£'K,iso,52 is in units of 10^^ erg and «o is in cm"-*. As 
noted in Section|231 if the marginal radio detection is indeed 
real, then this upper bound can be replaced with an equality. 
Finally, we can place a lower limit on the density using the 
condition that Vc < vx (i.e., v^ < 2.4 x 10"' Hz; 0.1 keV) 

no > 4.5 x 10-4(1 +z)-''hf^E2ils2 cm'^ (6) 

It is clear that ZiK.iso and no depend sensitively on our choice 
of z, £(. and eg- The fractions e^eg are not expected to exceed 
^ 1/3. We therefore calculate /iKiso for two cases: 1: = 
eg = 1/3, and 11: more typical values of e^, = = 0.1. We then 
calculate the range of allowed no set by Equations |5] and |6l 
which becomes wider as the redshift increases. For Case 1, 
this requires that z > 0.2, below which the constraints on no 
conflict (Figure |4|l. At the median observed redshift of the 
short GRB population, z ^ 0.5, we obtain /iK.iso « 3.7 x 10^" 
erg and no = 0.01-0.06 cm"^. For Case 11, the constraints on 
no require a higher redshift of z > 1 .5 (Figure|4]i. For a fiducial 
redshift of z = 1 .5, we obtain /iK.iso ~ 1 .2 x 10^^ erg and no = 
0.008 cm"^. The parameters for the two cases are listed in 
Table |4] Although we cannot distinguish between these two 
scenarios, both cases require low circumburst densities of n ^ 
0.01 cm-l 

3.4. Jet Opening Angle 

In the context of a jet break, we use the time of the break 
from the X-ray light curve (2.0 ±0.5 days) and the circum- 
burst density and energy estimates from the previous section 
to constrain 9j. The t ime of the break is a direct reflect ion of 
the jet opening angle dSari et al.lll999l:lFrail et ani2001l) . 
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e, = o.u^^(i^zry^E^!l,^^ (7) 

where tj^i is expressed in days. For our fiducial Case I (z = 
0.5, = £5=1 /3), £K.iso « 3.7 X 10^° erg and n « 0.01 -0.06 
cm"^ give 6j = 7-8°. This leads to a beaming correction on 
the energy of ft = [1 -cos(6';)] = 0.007-0.01, and therefore 
a true kinetic energy Ek = fbEK,mo « (3-4) x 10'*^ erg (Ta- 
ble IHi. To estimate the beaming-corrected 7-ray energy, we 
infer E^ iso from the ^wZ/f/BAT fluence and apply a bolomet- 
ric correction factor of 5 to roughly convert to a representative 
observed 7-ray energy range of ~ 10- 1000 keV. This factor 
is derived from short GRBs ob served by satellites w ith wider 
energy coverage (.Berger,2010t iMarguttiet all2012l) . We ob- 
tain ii^.iso = 2.1 X 10^" erg and therefore a true 7-ray energy 
of E^Ki 2 X lO'^'^erg. 

For Case II (z = 1.5,e^ = = 0.1), where no « 0.008 cm"^ 
and /iK.iso ~ 1 .2 X 10^^ erg, we obtain a smaller opening angle 
oi 6j K^3°. This leads to /}, w 1.4 x 10"-^ and hence, Ey w 
3 x 10"^** erg and Ek~2x 10'*'* erg. 

In both cases, the true 7-ray energy is few x lO'*** erg while 
the kinetic energy is an order of magnitude higher at z = 1.5 
than at z = 0.5. This results in a lower 7-ray conversion effi- 
ciency {rj^ =E~f/Etot) for Case II of 0.15 compared to 0.3-0.4 
for Case I (Table |4]i. The total energy even for Case II is 

10- 100 times lower that for long GRBs. 

3.5. Extinction 

We investigate the presence of extinction by comparing 
the unabsorbed X-ray flux and the optical non-detection at 
St = 17.7 hour. Since we do not know the exact location of 
the cooling frequency, we assume a maximum value i^cmux 
of 2.4 X 10'^ Hz (1 keV) and exti-apolate the X-ray flux 
to the optical band using the shallowest possible slope of 
(3 = -{p- l)/2 = -0.55 to obtain the lowest bound on the 
expected optical afterglow flux in the absence of extinction; 
any other assumption for the location of Vc < vx would 
result in a higher predicted optical flux density. For /:> = 2.1 
we estimate ivopt ~ 1.1 /iJy (' = 23.8 mag). Given that 
our observed 3(T upper limit is / > 24.4 mag, this implies 
a lower limit on the optical extinction in excess of the 
Galactic value of A, > 0.6 mag'". In the rest frame of 
the burst for a Milky Way extinction curve, this translates 
to ^?^^>. 0.6 mag for z = 0.5 and A^"**' > 0.2 at z = 1.5 
(ICardelh et al.lll989h . Using Galactic rel ations between Nh 
and Ay, A^H.intMy (1.1-2.2) x lO^' (iPredehl & SchmittI 
119951; IWatsonll201lb . we infer lower limits of A^H.int > 10^' 
cm"^ at z = 0.5 and A^Hant > 4.4 x 10^" cm"^ at z = 1.5, 
consistent with our value of 7.5 x 10^' cm"^ (z = 0) derived 
from the X-ray spectrum (Table [T}. However, an appreciable 
extinction is unexpected given the burst's location on the 
outskirts of its potential host galaxy. We note that the burst 
is located at Galactic coordinates (l,b)= (359. 3°, -19.4°) 
which is to ward the Galactic Bu lge on a steep gradient in the 
dust map ( Schlegel et al.r i998) and thus may be subject to 
substantial 30%) uncertainties in the Galactic extinction". 
Taking this uncertainty into account reduces the required 

We note that for p < 1.9, no host galaxy extinction is required 
" Using a high-resolution (SpwHM = 15") WISE 12^m map, we do not 
see strong evidence for any thin dust filaments at the location of the burst 
which would result in > 30% uncertainties in the Galactic Ay (D. Finkbeiner, 
private comm.) 



^host > 0.2-0.3 mag depending on the redshift of the burst. 



4. DISCUSSION 

4.1. Environment 

From our broad-band observations, we constrain the cir- 
cumburst density of GRB111020A to no ~ 0.01 cm"^ 
which is consistent with the low values inferred for a 
few previous short GRBs (Soderberg et al. 2006; Panaitescu] 
2006; Stratta_et a l. 2007; Perley et al. 2009b; Berger 2010] 
iFong et akl |201 ih . The inferred density fits well with the 
framework of NS-NS/NS-BH binary progenitor systems, 
which may be subject to substantial kicks from their host 
galaxies and are predicted to typically occur at densit ies of 
~ 10"^- 1 cm"^ (Perna & Belczvnski 2002; Belczvnski eLalJ 
2006). 

GRB 1 1 1020A has an offset of «0.80" from its most prob- 
able host galaxy (G3; Figure |2|i. For redshifts between 
z = 0.5- 1.5, this translates to a projected physical offset of 
« 5 - 7 kpc, which is comparable to the median value of 
~ 5 k pc for well-local i zed short GRBs wi th host associa- 
tions (iFonget al.ll20Tot IChurchet al.ll2011h . Although G3 
has the lowest probability of chance coincidence by an or- 
der of magnitude (Figure |3]l, we cannot rule out the possi- 
bility that G3 is a faint star The next most probable galax- 
ies, Gl and G2, are situated 2.8" (17-24 kpc) and 6.5" 
(40-56 kpc), respectively, from GRB 111020A (Figure EJ. 
If the burst originated from one of these galaxies, this would 
put GRB 111020A at the high end of the observed offset 
distribution, similar to the growing sub-class of apparently 
"hostless" short GR Bs which li kely occur > 30 kpc from 
their host galaxies ( Bergerll2010l) . All of these inferred off- 
sets are consistent with predicted offset distributions of NS- 
NS /NS-BH binaries originating in M ilky Way-typ e galaxies 
(Bloom et al. 1999; Fryer etal. 1999 HBelczynski e t al. 20061 
Salv aterra et al.l i20 10).' 

Most short GRB host galaxies with confirmed spectro- 
scopic redshifts hav e measured luminosities of « 0. 1 - IL* 
(iBerger et aLlfeOOTh . The apparent magnitude of G3 is / w 
24.3, which corresponds to Lb « 0.1 - IL* over z ~ 0.5- 
2.3 when compared to the luminosity function of galaxies 
at corresponding redshifts in the DEEP2 an d LBG surveys 
(IWiUmer et al.ll2006HReddv & Steidelll2009h . This is consis- 
tent with the redshift range inferred from the afterglow. 

We next investigate the nature of the dust and gas in the 
environment of GRB 111020A through an analysis of Ajj,"**' 
and A^H.int- We have shown that the burst requires dust extinc- 
tion of Ay°'*' > 0.2-0.6 mag, depending on the redshift of the 
burst and the uncertainty in Galactic extinction. We have also 
measured a neutral Hydrogen column density intrinsic to the 
burst environment of A^H.int = (7.5 ±2.0) x 10^' cm"^ at z = 0, 
which becomes higher for any other choice of z. High values 
of both dust extinction and X-ray absorption have been linked 
to "d ark" GRBs (e.g. Perley et all l2009al: ICampana et aP 
1201 lb which have optically sub-luminous afterglows com- 
pared to their X-ray or NIR counterpart s and can quantita- 
tively be classi fied by \Pox\ 5, [fel~0 .5 (Ivan der Horst et al] 
120091: see also iJakobsson et al.ll2004 . Wifli \Px\ = 1-0 and 
\l3ox\ < 0.46, GRB 111020A is consistent with this defini- 
tion of dark GRBs. While optical extinction intrinsic to long 
GRB environments is not uncommon and commensurate with 
their origin in dusty, star-forming regions, evidence for sub- 
stantial extinction has been reported for only one other short 
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burst, GRB 070724A, which required A^)?'' > 2 mag to ex- 
plain the supp ressio n of optical emi ssion relative to the NIR 
dBer^er et alTl2009l: iKocevski et ani2010l) . The location of 
GRB 070724A on the outskirts of its host galaxy, ~ 5 kpc 
from the center, suggested either an origin in a star-forming 
region or a progeni tor system which produced the dust itself 
dBerger et al .112009b . The potentially appreciable extinction 
and the location with respect to its putative host suggests that 
the same conclusions may be drawn for GRB 1 1 1020A. 

On the other hand, the relation between A^H.int and the dark- 
ness of a burst is less clear. A recent study of long dark 
GRBs shows them to have higher intrinsic column densi- 
ties than non-dark GRBs, which suggests that the darkness 
of a burst is largely due t o absorption by circumburst mate- 
rial ("Camp ana et al.ll2()llb . To investigate this relationship for 
GRB 1 1 1020A, we extract spectra and best-fitting A^H.int for all 
short GRBs with XRT-detected afterglows in the same manner 
as GRB 111020A (see Section |2. 2. Il l, over time ranges with 
no evidence for spectral evolution. There are 22 short bursts 
with sufficient X-ray counts to perform spectral analysis, 1 1 
of which have known redshifts (Table |5]l. We find a short 
GRB weighted average of A^H,int(z = 0)= (1.1 ± 0.14) x 10^' 
cm-2 (90% C.I.; Figure|5]l. In comparison, GRB 1 1 1020A has 
a high value of A^H.im = (7.5 ±2.0) x lO^' cm'^ at z = (Fig- 
ure |5]l. Taken at face value, it is surprising to find a large 
NH.mt for a substantial offset, and may suggest that the burst 
occurred in a relatively metal-rich environment. 

4.2. Beaming, Energetics, and Rates 

We uncover a break in the X-ray light curve of 
GRB 1 1 1020A at w 2 days, which we interpret as a jet break 
(Section [3.21 ). Depending on our values for z, Ce and e^, we 
infer an opening angle of « 3-8°. This is reminiscent of 
the first jet brea k discove ry in GRB 05 1221 A, with Oj -1° 
jSoderberg et al . 2006; Burrows et al.l |2006), and suggests 
that at least a fraction of these events are highly collimated. In 
addition, temporal b reaks at tj < fe w hours in GRBs 061201 
jStratta et al.1 120071) and 0905 1 jDe Pasquale et al.1 120T0I; 
iNicuesa Guelben zu et al."2012'). if interpreted as jet breaks, 
lead to 9j ^ 1° (Figure |6]l. However, these two cases resem- 
ble early breaks in long GRBs that are generally attributed to 
the cessation of energy injection, and not collimation. 

Although the remaining short GRB afterglow data is sparse, 
the lack of observed jet breaks in their X-ray light curves 
can be used to place lower limits on the opening angles. 
Indeed, Chandra observations of GRB 050724A out to 22 
days indicated Oj > 25° , consistent with a sp herical explo- 
sion (Griipe et al.l I2OO6I) . A recent study by ICoward et al.l 
(120121) analyzed the sample of short GRB Swift/XRT light 
curves up to August 2011 with monitoring > 1 day which 
included 6 additional events, and inferred Oj > 6 - 16°, as- 
suming «() = 1 cm"^ for all bursts. We revise this anal- 
ysis for 3 of the events with robust X-ray light curves 
(GRBs070714 B, 070724A, 07 1227; data analysis prescrip- 
tions from Marg utti et al.l201^ employing a more representa- 
tive no ~ 10"^ cm"^ (e.g. 'Sode rberg et al. '2006 and this work). 
We derive is-y.iso from the reported fluences, applying a bolo- 
metric correction when necessary to represent an energy range 
of ~ 10- 1000 keV, and infer more realistic lower limits of 
> 2-6° (Figure |6]l. These limits are indeed lower than the 
detected values for GRBs 05 1221 A and 1 1 1020 A, and there- 
fore do not add strong constraints on the distribution. We cau- 
tion that the sample presented here represents only the ^ 30% 



of the Swift short GRB population that have bright X-ray af- 
terglows and relatively slow flux decline rates; the remaining 
fraction do not have detectable X-ray afterglows or fade too 
quickly so constraints cannot be placed on their collimation. 

There are now two short GRBs with opening angle mea- 
surements, two with measurements based on more tentative 
early breaks, and an additional four with lower limits (Fig- 
ure |6]l. These early constraints create a distribution that may 
mimic the distribution for long GRBs, whi ch ranges from 
~2-20° with a median of 7° (Fi gure H Frail et al. 2001; 
iBerger et al.ll2003t iBloom et al.ll2003; Ghirianda et al. 20041 
Friedman & Bloom 2005; Racusin et al. 2009; Cenk o et al.l 
| 2010l .Filgas et al...2011; .Goldstein et aLi2011i; iCenko et al.l 
I2OI lb . More events are needed to assess the real differ- 
ences between the distributions. However, simulations of 
post-merger blac k hole accretion predict jets with 9j ~ 5-20° 
(lAlov et al.l2005 ; Rosswog 2005; Rezzolla et al. 201 1) t o sev- 
eral t ens of degrees (Ruffert & Janka , 1999b; .Rezzolla et al.l 
I2OI lb depending on the mechanism of energy extraction 
and Lorentz factor, so there are expectations on theoretical 
grounds that the short GRB distribution is wider 

The first major ramification of collimation is the correc- 
tion to the total energy release: the true energy is lower 
than the isotropic-equivalent value by the beaming factor, 
For GRB 111020A with an opening angle of « 3-8°, this 
correction factor is substantial, 0.001-0.01. Depending on 
the redshift, the beaming-corrected energy of GRB 1 1 1020A 
is oi (2-3) X 10^*^ erg (Table HJl which is an order of 
magnitude l ower than for GRB 05 1 22 1 A with » ( 1 - 2) x 
10'*'' erg (Soderberg etiD l2006l iBurrows et all 120^61) and 
GRB050724A with E^ « (0.4-4) x 10^° erg (iGrupe et al.l 
l2006b . The three remaining events with opening angle lower 
hmits, GRBs070714B, 070724A and 071227, have ranges 
of E-^ !v 10'^^ - 10^^ erg, where the upper bound is set by 
the isotropic-equivalent 7-ray energy in the ss 10- 1000 keV 
band. The small population of short GRBs with measured E^ 
therefore has a median value of E^ ^ 10"*^ erg, which is an or- 
der of magnitude below Swift long GRBs (iKocevski & Butlea 
I2OO8; Racusin et al. 2009 ) and 2 orders of magnitude beloy v 
the pre-5w//T population ( Frail et al.l2001l:lBlooin et al.l2003b . 
Again, this sample is incomplete because we can only mea- 
sure Ej for bursts with well-constrained opening angles. 

In a similar vein, we compare the beaming-corrected ki- 
netic energy and total energy [Ek, ^tot) of GRB 11 1020 A to 
the values for other short bursts. Because /iK.iso is more sen- 
sitive to our choices for z, and eb, we infer different values 
for Case I and II. For Case I, we infer Ek ^ (3-4) x 10*^ erg, 
£tot = E-^ + Ek ~ (5-6) X 10*^ erg, and r/^ « 0.3-0.4. For 
Case II, we calculate Efc ~2x 10"*^ erg, istot ~ 2 x 10"*^ erg 
and 77.^ « 0.15. GRB 05 1221 A had £g » 8 x 10"*^ erg and 
a total energy release of p a 2.5 x lO'*^ erg ( ISoderberg et al.l 
l2006t iBurrows et al.]l2006b while GRB 050724 had a total 
energy of 10^°- 10^* erg. With £,ot ~ (0.5-2) x 10'*'* erg, 
GRB 11 1020 A may be on the low end of the total energy 
distribution, but more events with beaming-corrected ener- 
gies are needed to better characterize the distribution for short 
GRBs. 

The true total energy release of short GRBs has strong im- 
plications on the energy extraction mechanism. Two primary 
mechanisms, the therm al energy release f r om vi? annihilation 
in a b aryonic outflow (iJaroszynskill 1 9931: iMochkovitch et all 
1 1993b and magnetohydrodynamic (MH D) processes in the 
black hole's accretion remnant (e.g. iBlandford & Znajekl 
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fT977t iRosswQg et alJl2003h . give different estimates for the 
expected energy release. Predictions for vv annihilation are 
largely dependent on the mass of the disc and efficiency to 
produce pairs. Simulations of an outflow due to vv annihi- 
lation suggest beaming-corrected total energy releases could 
reach 10'*''- 10'*^ erg (.Ruffert & Janka 1 999b a: Popham et aT] 
1 999[ lRosswo3l2005l: iBirkl et al.ll2067l: iLee & Ramirez-Ruiz 
2007h . Higher energy releases can be obtained from MHD 



processes, which can produce luminosities of > 10^^ erg s ' 
10^" erg s~' when corrected for beam ing; Pophar n et al.l 
119991; iRosswog et all 120031; iLee & Ramir ez-Ruiz' 20071) de- 
pending on the nature of the magnetic field amplifica- 
tion. While the true energy releases of GRBs051221A and 
050724 A suggest that MHD processes may be powering these 
events (Bergeret al. 2005; 'Grupe et al.ll2006l: 'Soderberg et al.l 
|2006; BuiTows et al. 2006), the'total energy of GRB in020A 
is consistent with predictions for both scenarios. 

The second major consequence of beaming is that the true 
event rate is higher than the observed rate by the inverse of the 
beaming factor (i.e., /?true = ft^Robs)- Thus, beaming provides 
essential information for understanding the relation to vari- 
ous progenitor systems and is of particular interest since the 
NS-NS/NS-BH merger rates, which are a critical input for es- 
timates of Advanced L IGO gravitational wave detections, are 
highly uncertain (e.g., Abadie et al.ll20ldt iMetzge r & Berger 
|2012). The current estima ted observed short GRB volumet- 
ric rate is ^ 10 Gpc~^ yr~' (iNakar et al.l l2006). The estimated 
NS-NS merger rate is much higher: ^ 200-3000 Gpc"^ yr"' 
(iKalogera et alJl2004l:lN^ar et alJl2006l) . 

The discrepancy in these rates can be explained if 
short GRBs have typ i cal ^ , ~ 8° (/^' - 100; see 
also Metzger & Bergeij |2012|) . The determination of 
GRB 111020A's opening angle of 3-8° (fj;^ = 100-730), 
along with the small but increasing sample of opening angle 
constraints for short GRBs, implies that at least a fraction of 
these events are significantly beamed and that the true rate of 
short GRBs is at least - 100-1000 Gpc"^ yr"'. With a few 
additional opening angle measurements, this value can be sig- 
nificantly improved. Other proposed progenitor models, e.g., 
WD- WD mergers or accretion-induced collapse of a WD/NS 
(lOin etal.lll998HLev an et al. 2006; Metzger et al. 2008) have 
estimated rates of < 1000 Gpc "^ yr"' and ^ 0.1 - 100 Gpc"^ 
yr~', respectively (Lee & R amirez-Ruid 120071: iDarbha et al.l 
1201 Ol) . so if a large fraction of short GRBs have opening an- 
gles of < 25°, these systems may not contribute significantly 
to the progenitor population. 

5. CONCLUSIONS AND FUTURE WORK 

We have presented observations of GRB 111020A, utiliz- 
ing extensive coverage in the X-rays with 5'w;/f/XRT, XMM 
and Chandra to uncover a temporal break, most naturally ex- 
plained as a jet break. Our limit on the radio afterglow from 
EVLA combined with the inference that v^ < vx leads to a 
robust range on the circumburst density of ^ 0.01-0.1 cm~-^. 
The jet break time of « 2 days leads to an opening angle of 
3-8°, depending on the redshift and equipartition fractions, 
which leads to beaming-corrected energies of E-^ « (2 - 3) x 
lO"*** erg, Ek ~ (0.3-2) x 10^'' erg and £tot « (0.5-2) x lO^"* 
erg. This result, along with the previous jet break constraints 
for GRBs051221A and 050724A suggests that there may be 
a sprea d in true energy release, ^ ^P^" erg for short 

GRBs (Bergeret al. '2005^. lOrupe et al.ll2006l; ISoderberg et al.l 
I2OO6; Burrows et al. 20061) . 

Furthermore, our optical observations provide a limit on the 



afterglow and enabled the discovery of a putative host galaxy 
with / « 24.3 mag. A comparison of the X-ray and optical data 
at (5f = 17.7 hours provides a lower limit on the host galaxy 
extinction of Ay°*' > 0.2-0.6 mag. This is consistent with 
the high intrinsic column density from X-ray absorption when 
compared to the mean for the short GRB population. 

GRB111020A demonstrates that rapid multi- wavelength 
follow-up is vital to our understanding of the basic proper- 
ties of short GRBs: the geometry, energetics, and circum- 
burst densities. In particular, the search for jet breaks on 
timescales of > few days is imperative for placing meaning- 
ful constraints on the opening angle distribution. Ideally, the 
detection of breaks in both optical and X-ray data leads to 
an unambiguous and tight constraint on the opening angle; 
however, optical afterglows are only detected in ^ 30% of 
Swift short GRBs, while X-ray afterglows have been detected 
in ~70%. Furthermore, optical afterglows are intrinsically 
faint and subject to host galaxy contamination, making long- 
term monitoring highly challenging. Therefore, the jet break 
search is optimized in the X-ray band where the burst is not 
subject to such contamination and the afterglow brightness is 
virtually independent of the typically low circumburst densi- 
ties. The X-rays also allow for a measurement of the kinetic 
energy of the outflow. Deep radio limits provide additional 
constraints on the circumburst density and energy. The EVLA 
upgrade is now enabling us to probe events with relatively low 
energy scales of ~ lO'*^ erg and densities of < 10"^ cm~^. 

The collimation of short GRBs will undoubtedly further 
our knowledge of their true energetics and rates. While the 
former provides information on the explosion and energy 
extraction mechanisms, the latter is crucial for understand- 
ing the relation to various progenitor systems (e.g., NS-NS 
mergers). Significant improvement on t he estimated short 
GRB observed rate of - 10 Gpc"^ yr"' (iNakaret al.ll2006l) 
will have a critical impact on estimates for coincident short 
GRB-gravitational wave detections in the era of Advanced 
LIGOA'IRGO ( Abadie et alJ l2010). The uncertainty in the 
observed short GRB rate is dominated by the uncertainty in 
the beaming fraction and with only a handful of short GRB 
opening angles measured to date, the discovery of even a few 
additional jet breaks in the coming years will enable signifi- 
cant progress. 
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TABLE 1 

GRB 1 1 1020A X-RAY Spectral Fit Parameters 



Telescope 


Detector 


St 
(ks) 


H.int 
(10-2 




C-stati,/d.o.f. 


Swift 
XMM 
Chandra 
Swi//+XMM 


XRT 
EPIC-PN 
ACIS-S 
XRT+EPIC-PN 


0.08-60 
61.4-76.8 
250.5-268.5 
see above 


1.0±0.3 
65+* 

4+2 -^ 

-OA 


2.2 ±0.5 
2.0 ±0.4 
1 

^•^-0.8 

2.0 ±0.3 


0.86/188 
1.0/256 
0.32/661 
0.94/446 



Note. — " These values assume a Galactic column density of A^H,gai — 6.9 X 10^" 
cm~^ jKalberla et alj2003) . using an XSPEC model of tbabs x ztbahs x pow at z = 0. 
' Uncertainties correspond to a 90% confidence level. 
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TABLE 2 
GRB 1 1 1020A X-RAY Observations 



5t 


Time Bin Duration 


Unabs. Flux (0.3-10 keV) 




(s) 


(erg cm ^ s ') 






Swift/XRT 




6.18 X lO'" 


7.44 X 10" 


(2.38 ±0.79) X 10- 


-10 


1.35 X 10^ 


4.64 X 10' 


(3.80 ± 1.03) X 10- 


-1 1 


2.66 X 10^ 


1.71 X 10^ 


(1.64 ±0.42) X 10- 


-11 


4.15 X 10^ 


1.26 X 10^ 


(2.43 ± 0.64) X 10- 


-11 


5.96 X 10^ 


2.36 X 10^ 


(1.03 ±0.26) X 10- 


-11 


7.97 X 10^ 


1.66 X 10^ 


(1.83 ±0.49) X 10- 


-11 


1 1 A •. , ^ a3 
1.14 X lU 


5.20 X 10^ 


(8.85 ± 1.72) X 10- 


-12 


5.94 X 10^ 


2.46 X 10' 


(1.65 ±0.33) X 10- 


-12 


1.17 X 10"^ 


2.46 X 10' 


(1.07 ±0.27) X 10- 


-12 


1.94 X 10"^ 


6.38 X 10^ 


(9.19 ±2.28) X 10- 


-13 


9 '^S V 1 n4 


6.28 X 10' 


(1.19±0.32)x 10" 


-12 


^ 1 Q Y 1 

j.i.y /V. lu 


5.90 X 10' 


(1.05±0.27)x 10- 


-12 


4.29 X lO'* 


1.61 X lO'' 


(8.36±2.41)x 10" 


-13 


1.26 X 105 


1.51 X 105 


(1.63±0.55)x 10- 


-13 


3.09 X 10' 


2.14 X 10' 


(l.ll±0.42)x 10" 


-13 




XMM/EPIC-PN 




6.91 X lO'* 


1.35 X lO'* 


(2.66±0.19)x 10" 


-13 




Chandra/ACIS-S 




2.61 X 10^ 


1.98 X lO'* 


(5.96±0.89)x 10" 


-14 


8.84 X 10^ 


1.98 X 10-* 


< 8.95 X 10-15 





Note. — Upper limits are 3o-. 
" These points were excluded from the broken power law fit. 
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TABLE 3 

GRB 1 1 1 020A Optical Photometry 



Date 


At 


Telescope 


Instrument 


Filter Exposures 


fpWHM 


Afterglow"'' 


p ab 


Gl" 


02° 


G3'' 


(UT) 


(d) 






(s) 


(arcsec) 


(AB mag) 


(/Tjy) 


(AB mag) 


(AB mag) 


(AB mag) 


2011 October 21.01 


0.74 


Magellan/Clay 


LDSS3 


r 3 X 360 


0.62 


> 23.4 


< 1.56 


> 23.4 


21.12±0.09 


> 23.4 


2011 October 21.01 


0.74 


Gemini-S 


GMOS 


i 9x180 


0.74 


> 24.4 


< 0.63 








2011 October 22.01 


1.74 


Gemini-S 


GMOS 


(• 11 X 180 


0.67 












2011 October 21.01+22.01 




Gemini-S 


GMOS 


i 20x180 


0.72 






23.89 ±0.17 


21.91 ±0.05 


24.27 ±0.16 



Note. — " These values have been corrected for Galactic extinction Ax iSchlafly & Finkbeiner|2Ql'T|) . 
Limits are 3cr. 
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TABLE 4 

Physical Parameters of GRB 1 1 1020A 



Parameter Case I [z = 0.5,ee = = 1/3] Case II [z = 1.5,ej = eg = 0.1] 



tj 


2.0 ±0.5 days" 


2.0 ±0.5 days" 


^7,iso 


2.1 X lO'^erg 


1.9 x 10^' erg 


^K,iso 


3.7 X lO'" erg 


1.2 X 10'^ erg 


no 


0.01-0.06 cm-3 


0.008 


Oj 


7-8° 


3° 


h 


0.007-0.01 


0.001 




2 X lO"** erg 


3 X lO'*^ erg 


Ek 


(3-4)xlO''8erg 


2 X lO"*' erg 


Etot 


(5-6)xlO''8erg 


2 X lO"" erg 




0.3-0.4 


0.15 


Note. 


— " Uncertainties correspond to a Icr 


confidence level. 
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TABLE 5 

Intrinsic X-ray Column Density of 

HYDROC.IiN. iVh im I'OR Smft SlIORT GRBS 



GRB 


z 


A'H.int 

(lO^i cm-2) 


a above zero 


050724 


0.258 


•'• -0.86 
< U.j4 


5.7 


05 1221 A 
060313 
060801 
061006 


0.547 

1.131 
0.438 


1 09+0.73 
n /I r+0 36 

0.45+;|;^* 
< 2.04 


4.5 
2.1 
2.4 


061201 




0.94+«|« 
< 1.89 


2.7 


0707 14B 
070724A 


0.923 
0.457 


4.2 


071227 
080123 


0.383 




6.8 
6.8 


080905A 
090510 


0.122 
0.903 


2.04+'j'|f 
< 0.80 


2.3 


090515 
090607 




<0.79 


3.2 


091109B 




< 1.58 




1001 17A 


0.915 


4.10+34' 
4.37!|| 
6.6l!|g 
<0.92 


2.2 


100702A 




2.1 


101219A 
110112A 


0.718 


3.3 



111020A ••• V-50!2^ 6.5 
111117A ••• 1-84!|^? 2.6 
llimA ■■■ 2.41+^?2 5.1 



Note. — Errors and upper limits quoted 
correspond to a 90% confidence level; z = is 
assumed when tlie redshift is not known. 
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Time after Burst 5t (s) 



Fig. 1 . — Unabsorbed X-ray flux light curve for GRB 1 1 1020A from Swift-XKT (red), XMM (green), and Chandra (blue). Flux errors are Icr. The Chandra So- 
upper limit is denoted by the blue triangle. The best-fit broken power law model (black solid line) for GRB 1 1 1020A is characterized by qj = —0.78, 02 = ~2. 1, 
and tj = 2.0 days. A single power law model with a = —0.78 (black dotted) violates the Chandra upper limit. Also plotted are X-ray light curves for short 
GRBs051221A (dark grey circles; Soderberg et al. 2006; Burrows et al. 2006) and 050724 (light grey circles; Grupe et al. 2006). The data for GRBs051221A 
and 050724 have been scaled for clarity. Grey lines trace the afterglow evolution with a break for GRB 051221A at 5 days and no break for GRB 050724A to 
^ 22 days. 
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Fig. 2. — Optical (-band observations obtained with GMOS on Gemini-Soutli. Left: Combined stack of two nights of GMOS /-band data. Stars Si and S2 
are labeled, as well as galaxies Gl and G2. X-ray positions of GRB 111020A are denoted by the circles (red: Swift-\KY, 1.6" radius, 90% containment; blue: 
Chandra, 0.33" radius, 90% confidence). Center. PSF-subtracted image with the centroids of SI and S2 (magenta circles). The subtraction reveals a third source, 
G3, with i « 24.3 mag. Right: Digital image subtraction of the two epochs obtained at 17.7 hours and 1.7 days, respectively, reveals no residuals in or around 
the Chandra position. 
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Distance from GRB 5R (arcsec) 



Fig. 3. — Probability of chance coincidence, P(< SR), as a function of angular distance from the center of the Chandra afterglow position for the three host 
galaxy candidates of GRB 111020A. The galaxy G3 has the lowest probability of chance coincidence P{< SR) = 0.01, and is therefore the most probable host of 
GRB111020A. 
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Fig. 4. — Isotropic-equivalent kinetic energy and circumburst density parameter space for GRB 111020A assuming = eg = 1/3 (left) and €e = eg = 0.1 
(right). The lower limit on the density (dark blue) is set by the condition that Uc < ux (Equation|6) while the upper limit (light blue) is set by radio observations 
(Equation|5). Also plotted are the values for Sk.iso at z = 0.1, 1 and 3 (red). The allowable parameter space set by these constraints is filled in yellow. 
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Redshift 

Fig. 5. — Excess neutral hydrogen column density, Wn int: versus redshift for GRB 1 1 1020A (1,2 and 3cr intervals denoted by green, red and orange lines) along 
with 6 Swift short GRBs with measured redshifts and optical afterglows (black filled circles) and 2 (GRBs 060801 and 101219A) with only X-ray afterglows 
(open circles). Also plotted are 1 1 short GRBs without secure redshifts (grey lines), 4 of which have only upper limits on A'H.ini (grey dashed). For GRBs without 
redshifts, the Nn i^t value at z = is scaled by (1 +z)^ * jGalama & Wiiersl2001h . En'ors and upper limits are at the 90% confidence level. The weighted mean for 
all short GRBs (black line) over the redshift interval z = 0— 1.2 is also shown. GRB 1 1 1020A has the highest WH.int of a short GRB to date and is well above the 
mean for short GRBs. 
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Fig. 6. — Distribution of opening angles for long (orange) and short (blue) GRBs. Aitows represent upper and lower li mits. The long GRB population includes 
pre-Svyf/? (Frail et al. 2001; Berger et al. 2003; Bloom et al. 2003; Ghirlanda et al. 2004; Friedman & Bloom 20(21), Swift jRacusin et al . 2009; Filgas et al. 201 1), 
and Fermi I Cenko et al. 2010; Goldstein et al. 201 1 ; Cenko et al. 201 1) bursts. The opening angle for GRB 1 1 1020A ranges from ~ 3 — 8° (depending on the 
redshift), while GRB 051221A has dj f» 7° iSoderbers et al. 2006; B urrows et ai.l2006l) . Tentative jet breaks (blue dashed) for GRBs 061201 (Stratta et al. 200^ 
and 090510 iDe Pasquale et al. 2010; Nicuesa Guelbenzu et al. 2012) are at ~ 1°. Shor t GRB lower limits are from the non-detection of jet breaks in SwiftfKSI 
data (this work, revised from lCoward et al.l2012l) and Chandra data for GRB 050724A IBerger et al.l200HGrupe et al.l200d) . 



